Molecular crystal approach for 7r-conjugated polymers: from PPP 
Hamiltonian to Holstein model for polaron states 



Stephane Pleutint* and Jean-Louis Fave* 
' Max- Planck- Institut fur Physik Komplexer Systeme, Nothnitzer Strafie 38, D-01187 Dresden 
* Groupe de Physique des Solides, 2 place Jussieu, 75251 Paris cedex 05 

(February 1, 2008) 

Abstract 

Starting from the 7r-electron Pariser-Parr-Pople (PPP) Hamiltonian which 
includes both strong electron-phonon and electron-electron interactions, we 
propose some strongly correlated wave functions of increasing quality for the 
ground state of conjugated polymers. These wavefunctions are built by com- 
bining different finite sets of local configurations extended at most over two 
nearest-neighbour monomers. With this picture, the doped case with one 
additional particle is expressed in terms of quasi-particle. Thus, the polaron 
formation problem goes back to the study of a Holstein like model. 
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I. INTRODUCTION 



The nature of the first excited states of conjugated polymers is an important and still 
unsolved question in condensed matter sciences p]. Knowing if they are band to band 
excitations or exciton states, if polarons, bipolarons or solitons are stable quasiparticles in 
the doping case, are fundamental issues for the understanding of the electronic properties 
of these compounds. 

The low-lying excited states are supposed to be suitably described by the well-known n- 
electron Pariser-Parr-Pople (PPP) Hamiltonian. This model Hamiltonian takes into account 
both strong electron-phonon and electron-electron interaction terms yielding only exact nu- 
merical solutions for the smallest oligomers 0. For the thermodynamic limit, the situation 
remains unclear since the calculations for the ground state and the excited states, including 
electron correlations, are uneasy to achieve and some drastic approximations are needed || . 

However, a first qualitative understanding of this complicated physics can be done by 
the use of some simplified Hamiltonian. For instance the Rice-Gartstein's molecular exciton 
model [|]], similar to previous works ||, is useful for a qualitative description of the linear 
absorption of conjugated polymers. On the other hand, the molecular Holstein's model gives 
a simplified picture of the polaron states . 

Recently, an approximate scheme to build the ground and the first excited states has 
been proposed ||. With this method, starting from the PPP Hamiltonian, one reaches a 
Rice-Gartstein like model; the excitations relevant for linear absorption are then easy to 
obtain and the results are comparable with those from more tedious methods ||. In this 
paper, we will show that the same procedure permits to derive formally, from the very 
same PPP Hamiltonian, the simple molecular Holstein's model for the polaron state in 
conjugated polymers. Polarons are thought to be important for the understanding of the 
charge transport in these compounds, and the possibility to study these non-linear states at 
a correlated level in an easy formalism is needed and valuable. 

We choose a simple dimerized linear chain as an effective model for conjugated polymers; 
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this chain is characterized by and r s , the double and single bond lengths, respectively. 
Extending our method to a realistic geometry is straightforward but the essential of the 
physics is reached with this simplified picture. Let us now briefly introduce the Pariser- 
Parr-Pople Hamiltonian which is our starting point 

Hppp = — ^ t n ^ n+ x{c^ „C n+ i >a + C^ +lj(7 C„ )(J ) + - ^2 Vn^C^Cn^ — -^){^ m>a iC m ,a' ~ ^) 
n i CT n,m,cr,a' ~ ~~ 

(1) 

where c\ la1 (c nj(J ) is the creation (destruction) operator of an electron in site n with spin 
a; t n ,n+i is the hopping term which includes via the electron-phonon interaction a linear 
dependence upon the length of the bond (n,n-hl) P,[10|. In the case of a dimerized linear 
chain, this dependence gives two distinct hopping integrals td and t s for the double and the 
single bonds respectively (| td \>\ t s |); they could be written as td/ s = t (l ± ^-5) where 
to is the hopping integral without dimerization, a is the electron-phonon interaction and S 
is a measure of the dimerization giving the difference of the lengths of single and double 
bonds [ID]. The Coulomb term is parametrized following Ohno, where the effect of the a 



electrons is considered via a phenomenologic screening, V nym = -, = == where r nym is 

y l+0,6117r n m 



the distance (in A) between two electrons localized on site n and m [IT|. We also write this 
term as V{r nym ) = V n ^ m and V = V{rd) for convenience. 

In view to link up the PPP Hamiltonian and the molecular crystal models, the Rice- 
Gardstein and Holstein models, we choose the monomer self-consistent orbitals as basis 
functions - this is the so-called exciton-basis ||. This choice is of course led by the dimer- 
ization. In our case, the monomers are the double-bonds and their self-consistent orbitals 
are associated with the following creation (destruction) operators for the bonding and anti- 
bonding orbitals: B® a = ^(4^ + c2+i jCT ) and A®, = ^(4^ - c2+i, CT ); here n indexes 
the double bonds. 

With this specific choice of local basis operators, the electronic configurations are built 
by combining different kinds of local configurations (LC) P,|T2"|]. In order to get a tractable 
model, we truncate the Hilbert space by choosing a small set of different LC which will 



be the elemental building blocks for the electronic configurations |L2]. These LC are the 
so-called generative local configurations (GLC) in ||. 

We may notice that this method shows some similarities with the Valence Bond method 



used efficiently for the studies of oligomers [|13|] but with the important difference that 
atomic sites are replaced by monomer units with internal electronic structure (double bonds 
here). The configurations build from GLC are diagonal with respect to the hopping term 
td on contrary to the Valence Bond configurations which are diagonal with respect to the 
Coulomb term. Each GLC is a set of several Valence Bond diagrams, chosen to be the 
adequate ones for a reasonable description of polymer states. 

In this work, we first improve the proposed ground state of ref. || by enlarging the set 
of electronic configurations used to describe it (section II). Second, we consider the case 
with an extra electron on the chain and show that, if one authorizes small lattice distortions 
around the extra particle, our treatment allows to reach quite naturally a Holstein like model 
but expressed in terms of many-body particle states (section III) . 



II. THE GROUND STATE 

We keep as GLC for the ground state the LC which appear the most relevant in calcu- 
lations performed on small oligomers ||. In ref ||, only three LC were considered; they 
are named F-LC, D-LC and Ct^-LC and are schematically represented in figure (|l|.a). This 
approximation could appear rather abrupt, but it is sufficient to get a correct qualitative 
picture of the linear absorption spectra as it was shown in ||; moreover, even at this level 
of approximation, the results are quantitatively comparable with the results of more te- 
dious calculations 0. In this work, we propose some natural improvements to this first 
approximation by extending the set of GLC. 

In a first improvement, we add to the previous set of GLC, the so-called Triplet-Triplet 
LC, TT-LC, shown in figure (|l].b), where two nearest-neighbour (n.n.) localized triplets are 
combined into a singlet. Together with the three first LC, they are the major constituents of 
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the ground state wave function in small cluster calculations 0. In a second improvement, 
we enlarge again the set of GLC by including in it the LC which interact directly with the 
four previous selected ones (figure|I[c). 

In the following, only the first case is treated explicitely. We develop in full detail our 
proposed way to get the ground state wave function with the four selected GLC. The case 
with the complete set of LC represented in figure ([!]) can be treated following the same 
scheme; only the obtained results are then given. 

First, we introduce the four GLC, their associated creation operators and their energies. 

• The named F-LC is associated with the creation operator 

= < T <! (2) 

This define the lowest LC in the range of parameter of interest; therefore we choose as 
reference state 

| >= n Fl | Vacuum > (3) 

n 

where | Vacuum > denotes the state without any ir electron. The state | > is the 
ground state considered in the molecular crystal approaches f|||; there, the linear 
dimerized chain is simply identified to an one-dimensional crystal of ethylene without 
any electronic correlations. 

With respect to | >, F\ = ft which is simply the identity operator. In the following, 
all the creation operators and the energies are defined with respect to | >. 

• The named D-LC is associated with the creation operator 

D\ L = A{^Ai A B^B nA (4) 

and with energy given by e<j = Atd- 
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The F and D-LC describe the dynamics of the 7r-electrons coupled by pairs into each 
monomer: the two electrons are independent in F-LC, whereas D-LC introduces in- 
tramonomer electronic correlation. In the strong dimerization limit, these two LC are 
sufficient to give a good approximation of the ground state; the system is then very 
close to a true molecular crystal. For small or intermediate dimerization, it is however 
necessary to consider more extended LC or, in other words, some fluctuations around 
the molecular crystal limit. This is done by introducing two more LC extended over 
two n.n. monomers. 

The named Ct^-LC is associated with the creation operator 

= l -(Al +1A B nA + Al +u B n>l - A^B n+u - Ai tl B n+u ) (5) 

and with energies given in the case of a linear dimerized chain by e ct = It^ + V — 
^{V(r s ) + 2V{r s + r^) + V{2rd + r s )). The last term, in bracket, is the attractive 
interaction between the electron and the hole due to the long-range part of the Ohno 
potential. 

The Ct5"-LC introduces n.n. intermonomer charge fluctuations, reproducing the con- 
jugation phenomenon in a minimal way. 

Last, the named TT-LC is associated with the creation operator 

^n,^ B nA A n+l,l B n+l,l + \,l B n,lA n+1 ^B n+1 ^ + A ] n ^B nA A ] n+lA B n+1A )) 

(6) 

and with energy given by e tf = 4t^ — (U — V) . 

In this LC, two Triplets appearing in n.n. monomers are combined into a singlet 
(figure([I].b)). It was shown to be important for the first time in the work of Schulten 
and Karplus (14j where it was recognized as a major constituent of the low-lying 
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excitations, the famous 2A g state, optically forbidden. In the ground state, which is 
our interest here, the importance of this LC can be comparable to the D-LC one 

We may notice that a similar treatment for the PPP Hamiltonian was proposed a few 
years ago to study the spin-charge separation mechanism in the limit of strong dimerization 

With our choice of four GLC, all possible electronic configurations are then build up. 
They are characterized by the number of D, Ctj~ and TT-LC, n d , n ct and n tt respectively, 
and by the positions of these different GLC. The positions of the D, Ctjf and TT-LC are 
labelled by the coordinates z{k) (k = l,..,n d ), y(j) (j = l,..,n ci ) and x{i) {i = l,..,n u ) 
respectively. The necessary non-overlapping condition between LC is supposed to be fulfilled 
all along the paper - the LC behave as hard core bosons. The electronic configurations are 
then expressed as 

ntt n ct n d 

| x(l),...,x(n tt ),y(l),...,y(n ct ),z(l),...,z(n d ) >= UR RTT^Ct^Dl^ I > ( 7 ) 

i=\j=\ k=i 

The GLC are all neutral local configurations, therefore the energy of (0) is independent of 
the relative positions between LC and entirely determined by the number of each GLC. 

E(n a , n ct , n d ) = n tt e tt + n ct e d + n d e d (8) 

At this point, we have to mention an incorrect statement in || where it is saying that the 
energy of the configurations made of F, D and Ctj~-LC depends on the relative positions of 
the Ct^-LC. This statement is actually wrong, however, this simplification goes in favor of 
our treatment (indeed, it was not possible to do calculations with this statement and finally 
the energy (|8|) was also adopted in ||). 

The way we choose to diagonalize the PPP Hamiltonian in the reduced Hilbert space 
spanned by the electronic configurations (0) follows from 0. First, we reorganize the con- 
figurations (0). We make linear combinations from the states with n d D-LC, n tt TT-LC 
localized at sites x{l), ...,x{n tt ) and n ct Ct^f-LC localized at sites y(l), ...,y(n ct ). Since we 
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are interested, at the end of the day, only by the lowest state in energy (the ground state), 
we can consider only the linear combinations of highest symmetry 

n d n t t n ct 

x(l), ...,x(n*)MV, -Mnctln, >= , E II D l(k) II II TT l(i) Ct l(j) I > 



Gn d 00)} fe=l 1=1 J=l 

(9) 

where the summation is carried over the £7N-2(n tt +n ct ) possible configurations. The energy 
of these combinations is still given by (§). 

The states (P) interact between them by the following term 

< x(l), 2/(1), y(n ct ),n d | if PPP | ac(l), x(n^), 1/(1), y(n ct ),n d + 1 >= 

^/(ra d + l)(iV - 2(ra tt + n ci ) - n d )^ 

The tri-diagonal matrix, where the diagonal part is given by (|8|) and the off-diagonal 
part by ([U]) can be divided into sub-matrices characterized by n ct localized Ct5~-LC and n a 
localized TT-LC but with a variable number of D-LC, n d (n d = 0, ...,2(n c t + ntt))] these 
sub-matrices can be separately diagonalized and it is easy to show that the resulting lowest 
states are given by the following expression 

| x(l), x{n tt ), 2/(1), y(n ct ) > c = a^^^l^y/c^^ 

| x(l), ...,x(n tt ),y(l), ...,y(n a ),n d > 



with energy expressed as 

E c (n tt , n ct ) = n tt e tt + n ct e ct + (N - 2(n tt + n ct ))e c (12) 

where 

e c = 2t d -±y/lGt 2 d + (U-V)* (13) 

The coefficients a and b of (Till) are given by a = , u ~ v = and a 2 + b 2 = 1. With these 
expressions, the double bonds free of Ctj~- and TT-LC are correlated independently. The 
upper-script, c, in ([n]) is for correlated. e c is called intramonomer correlation energy. 
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The next step toward the evaluation of the ground state is to retain, among all the 
states resulting from the previous sub-diagonalizations, only the lowest ones given by (PP. 
This approximation is well justified since the energy difference between these states and the 
corresponding lowest excited ones is given by the quantity ^16^ + (U — V) 2 which is rather 
high for usual parameters with a value around lOeV. We then reorganized the states ( |TT|) 
into collective excitations of highest symmetry 

| n tt ,n ct > c = [CZ; + Z; nct CZ +nc T" E I -,x(n tt ),y(l), ...,y(n ct ) > c (14) 

{m(i),y(j)} 

still associated with the energy (|T2|) and where the summation runs over the 
Cn~+n^ nct Cnu +nct possible configurations. The ground state is then expressed as a linear 
combination 

I GS >= x n tt ,n ct I n tt ,n ct > c (15) 

n t t,n c t 

where the coefficients X ntuHct are determined by solving the following secular equation 

I{n tt , n ct - l)^^ + (n tt e tt + n ct e ct - 2(n tt + n ct )e c - E)X nturict + I(n tt , n ct + l)X nu>rict+1 
[n t t(n c t + l)]^n tt ^-t s X ntt _^ nct+l + [n ct {n tt + \)\'2n ct Y t s X rHt+\,n ct -l = 

(16) 

where 

r( , I , 2(n tt + n ct ) - 1)(N - 2(n tt + n ct )) 2 

I{n tt ,n ct ) = \ (n ct + 1) a t s (17) 

V N -ntt- n ct 

The equation ([IB]) is not solvable with the interaction term (0). Next, and last, we do 
an approximation on the term I(nu,n c t) by assuming 

I(n tt ,n ct ) ~ ]j (n ct + 1)( N ^ 1 - n tt - n ct )V3a 2 t s (18) 

This is a very good approximation of (|TTD, if the number of GLC extended over two 
monomers, n 2 = n tt + n ct , is small Consequently, this treatment will be justified if 
in the final wave function, the most important configurations are the ones with a small 



value of U2'-, this is actually the case as it can be seen from the work of ref. || and as it 
appears, at the end of the day, in this study. 

With the last simplification, the problem is mapped onto (N — l)/3 independent three- 
level systems. One write 



x = Jc m ?- 1)/3) c n » +n ctv 



with < 



(19) 



yn t t,nct/yntt+l,n c t 1 
Vn,u,nct I ynu,n c t+l = C 

where E takes the integer part, 7 and ( are real constants to be determined. Inserting this 
definition in (|16|) and after some algebraic manipulations one finds that the problem goes 
back to calculate the lowest eigenvalue, e, of the following 3 by 3 matrix 



V3a 2 t s 



\ 



V3a 2 t s e ct -2e c *ft 



(20) 



/ 



The ground state energy is then simply divided into two different components 



Eos = Ne c + 



N 



(21) 



The first part is the intramonomer correlation energy defined by the first subdiagonalization; 
it is obtained by correlating independently the N double bonds. The second part is the 
intermonomer fluctuation energy defined by the second subdiagonalization; it is obtained 
by considering (N — l)/3 identical and independent effective three level systems defined by 
the matrix (]2"0D . Finally, the ground state wave function is clarified by the following two 
equations 

2a 2 e - e tt 



a 2 U 
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c 



(22) 



e y/3 e 

The resulting wave function contains, as the energy, two different kinds of components: the 
first ones localize electrons by pairs in the double bonds; the second ones introduce n.n. 
intermonomer fluctuations, charge fluctuations by means of Ct7-LC and spin fluctuation by 
means of TT-LC. 
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The ground state proposed above may be easily improved by adding new local configu- 
rations extended over two n.n. double bonds. For example, one can include the whole LC 
represented in the figure ([]]); the LC of (|l|.c) are the ones directly coupled to the others. The 
strategy is then the same. First, one takes care of the intramonomer correlation; second, one 
builds the collective excitations of highest symmetry; third, one approximates the part of the 
resulting interaction connecting configurations which differ by only one LC extended over 
two monomers in the way of ([181) . The problem is then equivalent to consider (N — l)/3 
independent seven-level systems; e is then the lowest eigenvalue of the associated 7 by 7 
matrix. 

In order to test our assumptions from which we propose several ground state wave func- 
tions in the form of (|T5|), we do comparisons, first, for the Su-Schrieffer-Heeger (SSH) model. 
For this model, similar to (|]) but without the complicated Coulomb term ||16|| , the exact 



result is well known pq , [17| . We compare this result with successively the results given by 
the model ground state of ref ||, the hereafter so-called model I, the one with in addition 
the TT-LC, the model II, and, last, the model with all the GLC represented in figure (|T]), 
the model III. We make comparison in function of the dimerization parameter x = ^-5. 
The results are shown in Table I where the percentage of the exact energy for our successive 
approximations are given. For x — 1, the case of complete dimerization, the three models 
give obviously the exact result. For x = 0, the case without dimerization, one gets around 
92% of the total energy. A-priori in this limit, one would expect less accurate results since 
the charge fluctuations of longer range than one play a role; they contribute in fact only in 
the missing 6%. For x = 0.15, a value often attributed to the polyacetylene, one gets around 
97% of the total energy. In conclusion, our approximation seems rather good for realistic 
cases, within this independent electron model. 

Next, we do also comparisons for the Hubbard model which is well known to be exactly 



solvable in one dimension this is the model (|I|) with a = and where only the on-site 
electron-electron interaction, U, is retained. For U = 0, one gets the SSH model without 
dimerization for which we obtained around 92% of the total energy (see Table I). Starting 
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from these values, the agreement monotonically decreases when U increases to finally get for 
infinite U, between 77% and 79% of the total energy, depending on the model (I, II or III) 
under consideration. This discrepancy shows that important LC are missing especially in 
the strong U limit; for instance, it is easy to see, just by energetic considerations, that, for 
large enough U, the TTT-LC, which is singlet made by three localized triplets, the TTTT- 
LC, which is singlet made by four localized triplets, and so on, may become important for 
the ground state wave function. With our specific choice of basis set completely localized on 
the double bonds, the dimerization parameter, x, is crucial; the more it will be important, 
the more our treatment will be relevant to become exact for a complete dimerization. In the 
Hubbard model, the dimerization is simply missing. If a ^ 0, the energy of the LC made 
from localized triplets increases making our approximations more and more reasonable. 

Last, we do comparison for the so-called extended Peierls-Hubbard model; this is the 
model (|l]) but with only the Hubbard term, U, and the n.n. interaction V, with the as- 
sumption that V = V"(rd) = V(r s ) pD| . On the contrary to the two previous models, 
this model is not integrable, also, we do comparisons with calculations performed with the 
Density-Matrix Renormalization Group (DMRG) technique P2 |; a very recent review of the 



advances related to this method may be found in |23||. The DMRG calculations have been 



done by E. Jeckelmann following the method developed in ref [20[|. We compare our 



approximate results with an extrapolation of the energy per unit cell made from calcula- 
tions for different lattice lengths up to two hundreds double bonds. The calculations are 
performed for a reasonable choice of parameters, U = 4i and V = t$. The results for 
several values of the dimerization parameter are confined in table II. We see that the errors 
are always less than 20% and are around 13 - 10% for realistic parameters. In our opin- 
ion, the agreements obtained here are satisfactory considering the relative simplicity of the 
wave functions proposed in this work. Moreover, with these approximate wave functions, 
some analytical insights are now possible which is very new in this range of parameters, 
appropriated for conjugated polymers. 

We do not compare for the moment our results with calculations made for the com- 

12 



plete PPP Hamiltonian. Nevertheless, since the remaining long range terms of the Coulomb 
potential are of smaller importance than the other terms of the Hamiltonian, one can rea- 
sonably expect only small quantitative changes in the results obtained with the extended 
Peierls-Hubbard model by using the full PPP Hamiltonian. 

Before closing this section, we may say that our wave functions are not variational since, 
in the way we choose to diagonalize the model, we do two successive sub-diagonalizations 
with some approximations. However, it is possible to build some variational wave-functions 



very similar to (|15). By the way, works are already done to propose a variational version 
of the model II |24] and other are in progress for the model III ||25|| . In other way, a very 
efficient Matrix-Product- Ansatz is also proposed in [p6| |. Compared to the work developed in 



24], one can say that our proposed way to diagonalize the PPP-Hamiltonian in the selected 



sub-Hilbert space is a very good approximation for appropriated parameters. 



III. POLARONIC STATES 

In this part, we consider the situation with one additional charge. We treat explicitely 
the case of an additional electron but, the case of the removal of one electron can be treated 
exactly in the same way. We show that this problem can be describe, with some approxi- 
mations, in terms of quasi-particles which obey to simple effective Hamiltonian. For a rigid 
lattice, we get a one dimensional tight binding Hamiltonian. If one authorizes some dis- 
tortions of the lattice around the extra particle, we get at second order in the distortion 
coordinates, a Holstein like model ||. In both cases, the parameters of these one-electron 
models are related to the PPP one's. 

In this work, we are not attempted to derive quantitative results. Our goal, based on 
semi-quantitative results, is to open up a way between a true many-body model given by the 
PPP-Hamiltonian and more simple one-electron models as the Holstein's model for polaronic 
states. Because it is not possible to solve the PPP model and since the important physical 
ingredients for an understanding of conjugated polymers are still not fully recognized [|TJ, 
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the derivation of more effective models is needed in order to get some physical insight. This 
work, and the very related one of ref. ||, goes in this direction. 

For convenience, we choose in this part the simplest description for the ground state 
given by the model I, using F, D and Ct^-LC. Since the model I already contains the most 
important local constituents for the ground state wave function, namely the F and CtJ -LC, 
we believe the results would not changed dramatically with a better description - by using 
the model II or III. Then, if we define 



| n d , n ct >= [C*--C d - 2n T 1/2 E Ct W- C i(n t) D\ {1) ■ ■ ■ D\ (nd) | > (23) 

{y(i)^O')} 

where the summation is over the C^~ nct C^~ 2rict possible configurations, the ground state 
wave function is simply written as 

N ct , N-2n ct 



\GS>= J2 <r nct K?y/C& E a N c - 2n ^b n c ^Cl~^\n d ,n ct > (24) 

n ct =0 n d =0 

where N a = E(^), a c = 4 + 6* = 1, a d = and «> + % = 1. e 



is then the lowest eigenvalue of the 2 by 2 matrix obtained from Q2"0D by suppressing the 
effective level corresponding to the TT-LC ||. For a typical choice of parameters relevant 
for conjugated polymers JIU| , the most probable LC is the F-LC (a 2 c ~ 0.98 and a^ t ~ 0.25); 



typical values for the energies are given by e c ~ —0.266^ and e ~ — 1.26eV. 

An additional charge disturbs the electronic cloud more or less strongly depending on the 
system under consideration. It could be a local distortion where the extra particle rearranges 
the system in short distances to create around it what it is called polarization cloud; this is 
the case for usual semi-conductors. On the contrary, it could be a complete rearrangement 
of the system as for strongly correlated systems In our case, the first behaviour is 

concerned and a quasi-particle picture is reached. 

We describe the perturbations caused by the extra-particle - the polarization cloud - by 
introducing a new set of LC more or less extended, which we call Charged Local Configu- 
rations (C-LC); the term "charged" means that they contain explicitely the extra-particle. 
Some example of C-LC, extended over one, two and three double bonds are shown in figure 
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(0) where the extra-electron is represented by the thick arrow. In the case of a "macroscopic" 
rearrangement of the electronic structure - as it could be the case for strongly correlated 
systems - the maximum extension of the relevant C-LC would be of the order of the system 
size. In our case, this critical size is of the order of some monomer units only. 

All around these C-LC, we assume the electronic structure unchanged with respect to the 
ground state; therefore, we consider such charged configurations - strictly speaking, these 
are linear combinations of electronic configurations but we adopt the proposed terminology 
for convenience - 

| a n >=\ N L > <g> | C« > ® | N R > 
I Pn,n+1 > = \ N L > ® | C£ n+1 >®\N R -1> (25) 

I ln,n+l,n+2 > = | N L > O | C^ n+1 , n+2 > ® | N R - 2 > 

where, | C" >, | Cf „ +1 > and | C^ n+1)n+2 > are some C-LC extended over one, two and 
three nearest-neighbour double bonds respectively, | Nl > (| N R >, | N R — 1 >, | N R — 2 >) 
is the part on the left (right) of the C-LC, described in the same way as | GS >. With this 
crude description, a C-LC acts as a dramatic boundary which simply interrupts the chain: 
the system is separated into two chains both described exactly as the ground state; the 
boundary contains explicitely the extra-particle within a defined C-LC. The more extended 



C-LC are inserted in the ground state in the same way as (|25l). 

With our approximation, the energy of each charged configuration (as (|25|)) is given by 
the addition of two different terms. The energy of the isolated C-LC and the energy of 
the external parts to the left and to the right of the C-LC. Since | Nl > and | iV^ > are 
neutral, the external parts don't interact via the Coulomb potential with the extra-particle. 
However, the configurations (^) must be improved for more quantitative results. Indeed, 
in a better description, because of the presence of the P-LC, the relative weight of the F, D 
and Ct5"-LC, controlled by the coefficients a c , b c , a ct and b ct should depend on their positions 
on the chain. Moreover, with an additional particle, the electron-hole symmetry is broken. 
All the LC used in | Nl > and | N R > are in the same sector of symmetry - the proper one 
for the building of the ground state. This is the case, for instance, of the Ct^f-LC where the 
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charge transfer on the right and on the left are of the same importance. In the presence of 
the P-LC these two charge transfers are no longer equivalent; the symmetry is broken and 
this implies a coulombic interaction between the P-LC and the external parts. These effects, 
not considering in this work, would certainly modify the polarisation cloud in a sensitive 
way. One can say, in other words, that the 'embedding" of the C-LC, due to the part | Nl > 
and | Nr > of (|25|) , are not treated efficiently in this work. We believe this is here the main 
point to be improved in the future for more quantitative results. 

For usual values of the PPP model, one kind of charged configurations is smaller in 
energy than the other and in such way that a perturbative treatment is possible to do. 
These configurations are due to the C-LC, referred as the P-LC hereafter (P stands for 
Particle), associated with the following creation operator 

J*, = Al„Fl (26) 

and represented in figure (@.a). The extra-particle is immersed in the reference vacuum and 
gives the following charged configurations 

| n >=| N L > ® | P n > <g> | N R > (27) 

where n referred to the position of the P-LC, and with an energy given by £ n = e n + (N r + 
N l )e c +(N r + N l -3)%, with N r + N l = jV-1 and e n = t d + f + the energy of the isolated 



P-LC. By comparing with (^T|), we see that there is a loss of intramonomer correlation energy 
and a loss of intermonomer fluctuation energy with respect to the ground state; indeed, the 
additional electron occupies a site in which one cannot place D and Ct7-LC. This loss of 
energy is more important for the more extended C-LC. 

In the following, we consider explicitely only the charged configurations ( p7| ) since the 
effects of the other charged configurations can be taken into account by perturbation. With 
our approximation, because of the n.n. hopping integral, the P-LC can hop on the lattice 
with the help of the F-LC or the Ct5~-LC. With the former, the P-LC can hop from site to 
site on the monomer lattice (see figure Q). 
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< n | H ppp | n ± 1 >= J = a^ctir (28) 

In this expression, the product a^ct gi ves the probability to find a F-LC in the wave function 
(HI); the factor 1/2 in (^8|) comes from our choice to work with the monomer orbitals. 
Moreover, with our approximation, there exists also a n.n.n. hopping process with the help 
of the more extended GLC, the CtJf-LC (see figure 

< n | H ppp \n±2>= a 2 Jj (29) 

The additional factor 2 in the denominator comes from the fact that only one term from 
the CtJ~-LC (see equation (§)) is involved during the transfer; the coefficient a 2 ct gives the 



probability to find a Ct5~-LC in the ground state wave function (24). The n.n.n. transfer 
is of course less important than the n.n. one's. With the values for the parameters we use 
here, the values of these two hopping processes differ by one order of magnitude. Therefore, 
we neglect the n.n.n. effective hopping term in this work. 

The extra-particle (P-LC) can be dressed by perturbation. Some effects of the other C- 
LC appear then in renormalized energy and n.n. hopping term for the extra-particle. This 
dressing of the P-LC can be simply done by a second order perturbative treatment, giving, 
in one hand, the so-called polarization energy 

e P = E J^V (30) 
and, in the other hand, some corrections for the n.n. hopping integral J 

Jeff = EwAj-zt + y^*) ( 31 ) 

g gi C"n &5 &n 

In these expressions tg and t s > are some interacting terms between the P-LC and other C- 
LC. The inequalities | — - |<< 1 are respected for the values of the parameters we use 
which guaranty the relevance of a perturbative treatment. After that, we have reached a 
quasi-particle picture, the quasi-particle being represented by the P-LC. 

In principle, many C-LC give some contributions to the perturbative series ( |5UD and (^TJ) . 
However, because the states (|25|) ignore many effects due to an inappropriate embedding, 
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as we already mentioned, we believe it is not useful to carry out the full calculation. Con- 
sequently, we do here a simplified treatment for the dressing of the extra-particle which we 
believe contains anyhow the most important contributions to (|30| ) and (|31|). This simplified 



treatment consists to consider the C-LC not embedded in the ground state defined by ( |24| ) 
but in a simplified vacuum made of only F-LC. Since, the F-LC is the very most important 
LC in the ground state (|24]) , we believe this simplified treatment sufficient to capture the 
most important parts of the polarisation energy and the effective hopping term. Moreover, 
among the remaining charged configurations only a few are incorporated in the perturbative 
treatment; they are shown in the figure (|2].b). By this last simplification we neglect all the 
C-LC shown in the figure (|5|) which take into account some long range polarisation effects; 
these C-LC are numerous but their total effect on fl30|) and (|31|) are small and they don't 
participate sensitively to the binding energy of the polaron state which is the main quantity 
we are looking for here. With this treatment, the corrections for the hopping term remain 
always negligible in the range of parameters of interest; we will therefore neglect these last 
corrections, J e ff- 

After the dressing operation, we obtain formally a one particle like problem with two 
characteristic energy terms E n , the site energy of the additional 'electron' with respect to 
the ground state, and J, the hopping term, which are functions of the PPP parameters: 
E n = e n — e c — |e + e p and J = ^. If we suppose a rigid lattice, the problem can obviously 
be diagonalized, giving a band centred at E n , with a bandwidth of 4 | J |. In the case 



of the SSH Hamiltonian [|T^] and by neglecting e and e p , the bottom of the band is given 
by E n =| td — t s |, the exact result; with inclusion of these corrective terms, this energy 
becomes slightly overestimated. The effective mass associated with the P-LC is given by 
m* ~ (a is the unit cell length) which is of course higher than the effective mass of a 
free particle on the bottom of the conduction band. With the Coulomb term, E n increases 
and m* stays unchanged. In conclusion, for a rigid lattice, we have reached a simple tight 
binding Hamiltonian - the so-called Hiickel model. Last, one may say that such approach is 
quite close in spirit of a recent work of J. Grafenstein et al., where an effective tight-binding 



model is derived at ab-initio level by means of an incremental method |27[ . 

Now we allow a relaxation of the lattice. For simplicity, we choose a model displacement 
where the two 'atoms' of the same double bond move with the same amplitude but in 
opposite directions (cf. figure (^[)). The two parameters E n and J depend now on the lattice 
coordinates via mainly the linear dependence of the two hopping terms, td(x n ) = td — ax n 
and t s (x n ,x n+ i) = t s + a(^f + ^y 1 )- On the contrary, the coulombic terms remain almost 
unchanged after a small displacement. The contributions due to these displacements to E n 
and J are small so we make a linear expansion with respect to {x n } of these two quantities 

E(x n ) = E n - a(a x n + ai(x n+1 + x n _i)) (32) 

J (x n ) = J - ab (x n + x n+x ) (33) 

where ao, a\ and bo are functions of PPP parameters and a is the electron-phonon interaction 
term ]TJ|. The extra elastic constraint of the dimerized chain due to the lattice relaxation 
in presence of an additional charge is expressed as 

E el = \^K eq [xl + C-^ + X -^f] (34) 

where K eq , the spring constant, is defined relatively to the dimerized equilibrium structure. 

The coefficients C n of the Holstein polaron wave function |J, | *ff p >= J2nC n {x n ) I n >, 
are determined by minimization of the corresponding total energy, ET({x n }), with respect to 
the lattice coordinate x n . At the second order in x n and taking into account that ~ 0.1 A 
in conjugated polymers [[J, we obtain the characteristic equations of the molecular Holstein's 
model 

[Fx n - 2 J - e]C n + JC n+1 + JC„_x = (35) 

x n k = F | C n | 2 (36) 
where the coefficients are expressed in function of the PPP parameters: F = (ao+2ai+4&o)a, 



k = 2K eq and J = By injecting ([5g) into (|35|) we obtain the non-linear Schrodinger 
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equation which gives the coefficients of the wave function; the second equation connects in 
a simple manner these coefficients and the lattice deformation. The analytical solution of 
these two equations in the continuum limit |J, valid for the "large" polaron case, gives the 
well known polaronic wave function C n = ^sech(7(n — no)) with Eb = -f^-, i] 2 = -y and 
7 2 = Yi the polaron state is localized around n , an undetermined quantity because of the 
translational invariance of the system. The associated binding energy of the polaron state 
is given by E p = 

We evaluate these quantities for several choices of parameters by the following sequence 
of calculations. First we optimize the dimerized geometry referring to a spring constant, 



K, relative to a hypothetical undimerized geometry fllql ; then, we evaluate K eg , calculating 
the second derivative of Et with respect to the dimerization coordinate at the geometrical 
equilibrium. Second, we solve the equations and (p6|). 



In the continuum version of the SSH Hamiltonian limit p8j , analytical expressions 



have been given. Our results always overestimate the reported values. For example, with 
t = 2.5eV, a = AleVX' 1 and K = 21eVA~ 2 , we get E p = 0.1 leV in place of OMAeV. In 
the same manner, our method also overestimate the value of the dimerization. These over- 
estimations occur naturally from our starting point which relies on a molecular description. 
Besides it has been shown that the SSH Hamiltonian is never equivalent to the Holstein's 



model for the dimerized linear chain ||28|| , so that the approximations of our model cannot 
be expected to lead to a good agreement in this case. However our approximations will 
cope better when the Coulomb interaction is taken into account; then the energies of charge 
fluctuation components decrease with their extensions, due to the long range part of the 
potential. This fact is in favour of our approximation. Furthermore the value of K used in 
this example is the appropriated one for the SSH Hamiltonian [[nj, but seems not to be in 
agreement with the experimental results obtained for small oligomers [[|. An higher value 
must be taken, favouring again our description. 

If one adds the Ohno potential, the binding energy decreases: as example, for the same 
choice of parameters and U = 11.16ey, we get E p = O.OOleV^. Finally, taking the same 
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parameters but with a more appropriate value for K, K = 41eVA~ 2 , we get a reasonable 
equilibrium geometry characterized by = 1.33A and r s = 1.47A Moreover we get the 
following values F ~ 9.5eVA~ 1 , J ~ l.leV, k ~ 78eVA~ 2 and the binding energy for the 
polaron decreases to E p ~ 0.025eV. In any case, we stay around traditionally adopted 
values. 

Before closing this section, note that with such a low binding energy, expected for con- 
jugated polymers, the quantum fluctuations of the lattice should be explicitely considered. 
However, it is for the moment totally hopeless to introduce additional bosonic variables in 
the full PPP Hamiltonian. 

IV. CONCLUSION 

In conclusion, we have proposed a simplified treatment of the PPP Hamiltonian which 
is typically a diagonalization of this Hamiltonian in a restricted Hilbert space. The method 
adopted, using monomer orbitals, is a natural way to bridge the gap between small cluster 
and polymer calculations |2|f|. The ground state is composed of intermonomer nearest- 
neighbour fluctuation components in the background of coupled electrons by pairs localized 
on monomers. Comparisons with DMRG results for the extended Peierls-Hubbard model 
show satisfactory agreements considering the simplicity of our proposed wave functions. The 
electronic excitations are then described as local perturbations moving in this "vacuum". 
For an appropriate set of parameters, this description gives rather good values for the dimeri- 
sation and for the energies of the excited states active in one photon spectroscopy |§. In 
the doped case (2N+1 particles) studied here, following the adiabatic scheme proposed by 
Holstein ||, we show that our model leads naturally to a Holstein's polaron like problem. 
However our description differs drastically from the Holstein's polaron image in the sense 
that it is able to describe the behaviour of a strong correlated (2N+1) particle state whereas 
the Holstein's model considers only the additional particle in interaction with a deformable 
medium. The obtained binding energy of the polaron is of the correct order of magnitude. 
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Some improvements are suitable concerning, first, the ground state where more extended 
GLC must be considered in order to reproduce more accurately the derealization proper to 
7r-systems. In other hand, variational calculations based on the very same ideas are possible 
24|^6[. For the doped case, we believe the first point would be to improve the description 



of the vacuum in presence of the extra-particle. Even if it is difficult to test our derivation 
in part owing to a lack of accurate calculations including correlation effects, we think that 
our formulation keeps the essential behaviour of the considered physical phenomenon and 
believe that it could be useful for future more advanced studies in part because of its relative 
conceptual simplicity and of its ability to give analytical expressions. For example, the 
behaviour of polaron states in presence of a strong electric field 0, which corresponds to a 
common situation in electroluminescence studies, would be to consider taking into account 
the effects of the strongly correlated N particle system. 
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FIGURES 



a) 



b) 



TT 



c) 



FIG. 1. Generative-Local-Configurations (GLC) selected to build the ground state wave func- 
tions. The set of GLC (a) defines the model I, the set of GLC (a+b) defines the model II and the 
whole set of GLC (a+b+c) defines the model III. 



a) 



b) 



c) -4- U + 

FIG. 2. Examples of Charged-Local-Configurations (C-LC) extended over one (a), two (b) and 
three (c) double bonds. The extra-particle is represented by the thick arrow. The quasi-particle is 
identified with the P-LC (a), the lowest C-LC in energy. 
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FIG. 3. Nearest-neighbour hopping process for the P-LC assisted by the F-LC 



t 



t 



FIG. 4. Next-nearest-neighbour hopping process for the P-LC assisted by the Ct 1 -LC. 




FIG. 5. Example of Local Configurations including long range polarization effects 
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n-1 n n+1 

FIG. 6. Model for the lattice deformation adopted in this work. 
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TABLES 

TABLE I. Percentage of the exact energy obtained with the different models studied here for 
the S.S.H. Hamiltonian. The model I contains the F, D and Ctj~-LC; the TT-LC are added for the 
model II and the whole set of LC shown in figure (|l|) are considered for the model III. 

model I model II model III 

x = 0. 91.6% 92.1% 92.7% 

a; = 0.15 96.7% 96.9% 97.2% 



TABLE II. Energy per unit cell for an infinite lattice obtained with the three successive approx- 
imations (model I, II) HI) an d DMRG calculations for the extended Peierls-Hubbard model with 
U = At and V = t; in the case of DMRG, the energies per unit cell are obtained from extrapolation 
of large cluster calculations up to 400 sites. 



x = 0.05 
x = 0.15 
x = 0.25 
x = 0.75 



model I 
0.373311 
0.306566 
0.236707 
-0.160370 



model II 
0.369217 
0.304317 
0.235678 
-0.159835 



model III 
0.366820 
0.303046 
0.235022 
-0.159840 



DMRG 
0.313599 
0.270381 
0.213969 
-0.164925 
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